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ABSTRACT 


A  statistical  approach  for  estimating  intervals  of  certification  or  biases  of  facilities 
or  measurement  systems  including  uncertainties  is  set  forth  based  on  MxN-order  level 
testing,  which  is  defined  as  M  repetitions  of  the  same  N-order  level  experiment  in  M 
different  facilities  or  in  the  same  facility  with  M  different  measurement  systems.  In 
absence  of  reference  values,  mean  facility  or  measurement  system  used  for  assessing 
intervals  of  certification  or  biases.  Certification  or  biases  of  facilities  or  measurement 
systems  are  defined  as  processes  for  assessing  probabilistic  confidence  intervals  for 
facilities  or  measurement  systems  for  specific  tests,  data  reduction  equations,  conditions, 
procedures,  and  uncertainty  analysis.  Similarly,  subgroup  analysis  performed  for 
isolating  and  assessing  levels  of  differences  due  to  use  of  different  model  sizes  (scale 
effects)  or  measurement  systems.  An  example  provided  for  towing  tank  facilities  for 
resistance  tests  using  standard  uncertainty  analysis  procedures  based  on  an  international 
collaboration  between  three  facilities.  Although  number  of  facilities  minimum,  the 
results  demonstrate  usefulness  of  approach  and  support  recommendation  of  future 
collaborations  between  more  facilities.  Knowledge  of  intervals  of  certification  or  biases 
is  important  for  design,  accrediting  facilities  or  measurement  systems,  and  CFD 
validation. 
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STATISTICAL  APPROACH  FOR  ESTIMATING  INTERVALS  OF 
CERTIFICATION  OR  BIASES  OF  FACILITIES  OR  MEASUREMENT 
SYSTEMS  INCLUDING  UNCERTAINTIES 

1.  INTRODUCTION 

Experimental  fluid  dynamics  (EFD)  testing  in  large-scale  facilities  at  research 
institutes  is  undergoing  change  from  routine  tests  for  global  variables  to  detailed  tests  for 
local  variables  for  model  development  and  computational  fluid  dynamics  (CFD) 
validation,  as  design  methodology  changes  from  model  testing  and  theory  to  simulation- 
based  design.  Detailed  testing  requires  facilities  utilize  advanced  modern  measurement 
systems  (MS)  with  complete  documentation  of  test  conditions,  procedures,  and 
uncertainty  analysis.  The  requirements  for  intervals  of  uncertainties  are  even  more 
stringent  than  required  previously  since  they  are  a  limiting  factor  in  establishing  intervals 
of  CFD  simulation  validation  [1]  and  code  certification  [2]  and  ultimately  credibility  of 
simulation  technology.  In  addition,  routine  test  data  more  likely  utilized  in  house, 
whereas  detailed  test  data  is  more  likely  utilized  internationally,  which  additionally 
requires  use  of  standard  procedures  and  uncertainty  analysis  and  establishment  of 
benchmark  intervals  of  uncertainties.  Detailed  testing  offers  new  opportunities  for 
research  institutes,  as  the  amount  and  complexity  of  testing  is  increased. 

Methodology  and  procedures  for  estimating  EFD  uncertainties  have  developed 
and  progressed  over  the  past  fifty  years.  Formalization  [3-5]  followed  by  standard 
procedures  with  emphasis  on  simplification  and  practical  application  [6-7].  However, 
rigorous  use  continues  to  be  a  problem  in  both  research  and  design  at  university,  industry, 
and  government  laboratories.  Another  problem  is  lack  of  methodology  and  procedures 
for  estimating  intervals  of  certification  or  biases  of  facilities  or  MS,  i.e.,  establishing 
intervals  of  confidence  for  facilities  or  MS  arising  from  systematic  errors  due  to 
differences  or  peculiarities  in  individual  facilities  or  MS.  Such  differences  or  peculiarities 
arise  from  detailed  facility  geometry  or  MS  design,  working  fluid  and  flow  quality, 
conditions  and  procedures,  test  engineers,  specific  locations,  etc.  Adding  to  the  problem 
is  fact  that  such  developments  require  considerable  resources  and,  in  case  of  facilities, 
cooperation  amongst  institutes,  which  often  crosses  international  boundaries.  Estimating 
intervals  of  certification  or  biases  of  facilities  or  measurement  systems  are  required  for 
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establishing  standard  intervals  of  uncertainties  for  various  types  of  facilities  (towing 
tanks,  wind  tunnels,  flumes,  etc.)  and  tests  (forces  and  moments,  motions,  waves 
elevations,  mean  velocities,  turbulence,  etc.)  and  MS  (load  cells,  potentiometers,  wave 
probes,  pitot,  LDV,  PIV,  etc.).  This  is  important  for  design,  accrediting  facilities  or 
measurement  systems,  and  CFD  validation. 

Most  work  on  facility  or  MS  biases  is  for  small-scale  flow  meter  calibration 
facilities  with  focus  on  validation  of  accuracy,  comparison  of  international  flow 
standards,  and  establishing  domestic  flow  traceability  [8].  Proficiency  testing  programs 
are  used  to  establish  flow  measurement  traceability,  which  are  largely  based  on  Youden 
plots  [9]  requiring  two  (e.g.,  tandem  and/or  upstream  and  downstream)  MS  at  each 
facility.  This  approach  not  easily  extended  to  large-scale  multi-purpose  facilities  with 
complex  MS,  including  consideration  individual  facility  and  measurement  systems  bias 
and  precision  limits.  Individual  facility  and  measurement  systems  bias  and  precision 
limits  are  required  for  use  of  such  data  as  well  as  helpful  in  MS  improvements. 

For  large-scale  facilities  such  as  wind  tunnels  and  towing  tanks  with  complex  MS, 
only  limited  work  done  and  facility  or  MS  biases  not  yet  considered.  The  NATO, 
AGARD,  Propulsion  and  Energetics  Panel,  Uniform  Engine  Testing  Program,  was  a 
remarkable  early  exercise  in  large-scale  testing  in  which  the  same  jet  engines  were  tested 
in  a  number  of  jet  engine  test  stands  in  various  NATO  countries  and  uncertainties  were 
estimated  to  explain  whether  data  scatter  was  within  the  data  uncertainty  and  conclusions 
were  drawn  [10].  Ref.  [11]  compares  results  from  wind  tunnel  tests  for  same  geometry 
and  conditions  at  two  different  institutes,  model  scales,  and  using  a  number  of  different 
measurement  techniques  and  extensive  error- analysis.  The  Cooperative  Experimental 
Program  of  the  Resistance  Committees  of  the  17-19  International  Towing  Tank 
Conferences  (ITTC)  [12]  compare  results  from  towing  tank  tests  at  22  institutes. 
Comparisons  are  made  of  global  (resistance,  sinkage  and  trim,  wave  profile,  wave  cut, 
wake  survey,  form  factor,  and  blockage)  and  local  (surface  pressure  and  boundary  layer 
traverses)  data  for  a  standard  geometry  (Series  60)  of  different  sizes  (1.2-9. 6  m). 
However,  uncertainty  assessment  not  considered.  The  cooperative  uncertainty 
assessment  example  for  resistance  test  of  the  Resistance  Committee  of  the  22nd  ITTC 
[13]  compare  results  from  towing  tank  tests  at  7  institutes  of  resistance  test  bias  and 
precision  limits  and  total  uncertainties  following  standard  uncertainty  assessment 
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procedures,  but  for  different  model  geometries  and  sizes  (Series  60,  container  ships,  and 
5415). 

In  the  following,  a  statistical  approach  for  estimating  intervals  of  certification  or 
biases  of  facilities  or  MS  including  uncertainties  is  set  forth.  N-order  level  testing 
reviewed  followed  by  definitions  for  MxN-order  level  testing,  which  defined  as  M 
repetitions  of  the  same  N-order  level  experiment  in  M  different  facilities  or  in  the  same 
facility  with  M  different  measurement  systems.  If  reference  values  known,  present 
approach  used  at  either  the  N-order  or  MxN-order  levels.  However,  unlike  CFD  where 
EFD  provides  reference  values,  for  EFD  reference  values  are  seldom  known,  e.g.,  from  a 
standard  facility  or  MS.  In  absence  of  reference  values,  mean  facility  or  MS  used  for 
assessing  intervals  of  certification  or  biases.  Herein,  certification  or  biases  of  facilities  or 
measurement  systems  are  defined  as  processes  for  assessing  probabilistic  confidence 
intervals  for  facilities  or  measurement  systems  for  specific  tests,  data  reduction  equations, 
conditions,  procedures,  and  uncertainty  analysis.  Similarly,  subgroup  analysis  performed 
for  isolating  and  assessing  levels  of  differences  due  to  use  of  different  model  sizes  (scale 
effects)  or  measurement  systems.  An  example  provided  for  towing  tank  facilities  for 
resistance  tests  using  standard  uncertainty  analysis  procedures  based  on  an  international 
collaboration  between  three  facilities. 

2.  ESTIMATING  INTERVALS  OF  CERTIFICATION  OR  BIASES  OF 

FACILITIES  OR  MS 

Designing  tests  for  estimating  intervals  of  certification  or  biases  of  facilities  or 
MS  requires  special  care  and  consideration.  Many  factors  affect  certification  or  biases  of 
facilities  or  MS;  therefore,  as  with  estimating  precision  limits  (i.e.,  random  errors),  only 
those  factors  specifically  isolated  (i.e.,  turned  on)  are  included.  For  example,  if  interest  is 
for  certain  types  of  measurements  using  same  MS  in  different  facilities,  then  model- 
geometry,  tests,  data  reduction  equations,  conditions,  procedures,  and  uncertainty 
analysis  should  all  be  the  same.  Ideally,  standard  models  are  used.  Otherwise,  effects  of 
differences  in  model  geometry  are  included.  If  models  are  geosyms,  but  two  different 
scales,  then  scale  effects  are  included.  Similarly,  if  different  MS  used  at  different 
facilities,  than  effects  of  MS  are  included  and  so  on.  Although  approach  used  for  either 
facilities  or  MS,  presentation  that  follows  is  for  facilities  since  same  as  for  the  example. 
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2.1.  N-order  Level  Testing 

In  N-order  level  testing,  N  repetitions  of  the  same  experiment  in  the  same  facility 
conducted 


1  N 

Xi=—Yx/ 

1  N  £  ' 


(1) 


Where  Xj  and  Xt  are  single  realization  and  individual  facility  mean  results,  respectively. 
The  uncertainty  in  Xt  is  given  by  the  root-sum-square  (RSS)  of  bias  B  v  and 
precision limits 


Ux=JBl+Pi 


(2) 


The  bias  limit  obtained  by  considering  all  sources  for  systematic  errors,  based  on  0-  or  1- 
order  testing.  The  precision  limit  obtained  by  the  standard  deviation  of  the  mean 


(3) 


Where  SW  is  the  standard  deviation  of  the  sample  population  Xj 


S V/  = 


1  N 

— Y(x/ 
N-\jr 


-xt)2 


1/2 


(4) 


Under  the  assumption  of  a  nonnal  distribution  for  the  sample  population  Xj ,  95% 
confidence  level,  and  N>10,  the  estimated  true  result  of  the  experiment  X ET^  lies  inside 
the  intervals 


and 


X!-L  <.V,,  <X/+U 


Xl 


(5) 


xi-uXi<xETN<xi+uXi 


(6) 
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for  the  single  realization  and  mean  experimental  result,  respectively,  where 

U\, =Bl+(2Sx,}  (7) 

Additionally,  this  assumes  facility  biases  J3FB  =0.0,  otherwise  XET  is  the  biased  estimated 

true  value.  XET  is  referred  to  as  estimated  true  result  of  the  experiment;  since, 

confidence  in  equations  (5)  and  (6)  relies  on  confidence  in  equations  (2)  and  (7).  Note 
that  at  the  N-order  level  outliers  often  discarded  if 

^/\  =  \x/-Xt\>2Sx{  (8) 


2.2.  MxN-order  Level  Testing 

In  MxN-order  level  testing,  M  repetitions  of  the  same  N-order  experiment  in  M 
different  facilities  conducted 


_  i  M 

X  = 


M 


M  1  M  N 

2A 


i=\ 


Mxn 


(9) 


where  X  is  the  mean  facility  result.  The  uncertainty  in  A  is 

u-x  =  P\ + pi 


(10) 


The  bias  limit  of  the  mean  5-  is  the  average  RSS  of  the  M  bias  limits  B , 


l  M 


M 


(11) 


;=i 


The  precision  limit  is  the  standard  deviation  of  the  M  results  Xt 


=^Sx  =  2—22=  =  —= 

4m  4m 


M- 


M  _ 

7Z«  --o2 
^  1=1 


ni/2 


(12) 


or.Fr,  is  the  average  RSS  of  the  M  precision  limits  Pv  from  equation  (3) 


/A  = 


i  M 

-L  yp2 

x' 


(13) 
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Under  the  assumption  of  normal  distribution  for  the  sample  population  Xt,  95 % 
confidence  level,  and  M>10,  the  estimated  true  result  of  the  experiment  XETm  lies  inside 
the  intervals 


x.-ux<xFT  <x.+ux 

l  Aj  MN  1  Xj 


(14) 


and 


X-Ux<XFT  <X  +  Ux  (15) 

a  ^LMN  A  v  7 

for  the  individual  and  mean  facilities,  respectively,  where 

Uli=Bl+(2Sx^)  (16) 

Equation  (16)  equals  equation  (2)  if, S'.  =  S  /  -Jn  .  Additionally,  this  assumes  the  mean 

facility  biases  ftFB=0.0,  otherwise  XET  is  the  biased  estimated  true  value.  Fig.  1 

displays  MxN-order  level  testing,  including  individual  and  mean  facility  results  and  their 
bias  and  precision  limits  and  total  uncertainties;  biased  parent-population  mean  value  p; 
and  estimated  true  experimental  result  XET  . 


2.3.  Using  Mean  Values  As  Reference  Values 

Di  defined  as  the  difference  between  the  N-order  level  individual  facility  Xj  and 
MxN-order  level  mean  facility  X  values 

Dt=Xt-X  (17) 


and  its  uncertainty  U D  is  defined  as  the  RSS  of  the  uncertainties  of  Xt  and  X 


uDl  =  Jul  +u| 


,  YbI  4  S 
U2D=Bl+Pl+^^  + 


M 


M 


(18) 

(19) 


6 


or 


ui=B‘+K  + 


Z<, 

M 2 


-  + 


10 

M2 


(20) 


If  the  absolute  value  of  Dt  is  less  than  U D 

\D,\  <UDi  (21) 

then  the  individual  facility  is  certified  at  interval  UD  ,  whereas  if  the  absolute  value  of  D. 
is  greater  than  UD 

\D,\>UDi  (22) 

then  the  facility  bias  UFB  is  defined  as 


with  total  uncertainty 


If  Dt  is  much  greater  than  U D 


U1  =D 2 

W  FB: 


U\  =U2Xi  +U 


2 

FB, 


»uDt 


(23) 


(24) 


(25) 


then  Dj  approximately  equals  the  individual  facility  bias  error  J3FB 


A  -  PFBt  (26) 

such  that  the  biases  can  be  estimated  in  both  sign  and  magnitude  and  used  for  calibration. 
Fig.  2  displays  use  of  mean  values  as  reference  values,  including  both  situations  of 
estimation  of  interval  of  facility  certification  and  bias. 

For  certified  facilities,  interval  certification  provides  additional  confidence  in 
accuracy  measurements;  since,  validates  X(  and  accounts  for  U-  in  assessing  level  of 

certification.  For  non-certified  facilities,  accounting  for  facility  biases  provides  UT  , 
which  is  an  improved  estimate  than  O' v  .  Presumably,  design  sets  the  requirements  on 


7 


appropriate  intervals  for  certification  or  biases  of  facilities.  Comparison  of  U x  with  U- 
and  uncertainties  from  other  facilities  is  useful  in  developing  strategies  for  reduction U x  . 
Note  that  for  sufficiently  large  M,  Ux  »U-  so U D  ~UX  .  In  this  case,  MxN-order 

level  testing  primarily  provides  X  .  As  already  mentioned,  if  reference  values  are  known, 
present  approach  used  at  either  the  N-order  or  MxN-order  levels,  which  for  completeness 
is  included  as  Appendix  A. 


2.4.  Subgroup  Analysis 

Isolating  and  assessing  levels  of  differences  due  to  use  of  different  model  sizes  (scale 
effects)  or  measurement  systems  is  of  importance.  Subgroup  differences  assessed  by 
comparison  of  the  subgroup  mean  to  the  total  mean,  with  consideration  to  the  uncertainty 
in  the  comparison.  For  L  subgroup  facilities,  subgroup  mean  and  uncertainty  given  by 


where 


(27) 


U-  =Bl  +P-2 

X  SG  X  SG  X  SG 


(28) 


(29) 


and 


p~x-=7ZiT^l;{Xs<i~Xsafr  (30) 

The  subgroup  difference  and  difference  uncertainty  given  by 


Dsg  =  X  sg  —  X 

(31) 

uL=uL+uj- 

(32) 

lAfG  -  UDse 

(33) 

differences  can  not  be  discerned,  i.e.,  within  noise  intervals  of  comparison,  whereas  for 
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(34) 


differences  discernable,  which  suggests  need  for  separate  certification.  Present  approach 
differs  from  analysis  of  the  means  [14];  since,  takes  into  account  both  bias  and  precision 
uncertainties  for  both  subgroup  and  mean  in  assessing  intervals  of  subgroup  differences. 
Analysis  of  the  means  compares  DSCJ  |  withf^  j  VZ  . 


3.  EXAMPLE  FOR  TOWING  TANK  FACILITIES 

It  is  not  easy  to  provide  an  example  of  the  proposed  approach;  since,  as  already 
mentioned  it  requires  considerable  resources  and  cooperation  amongst  institutes  often 
crossing  international  boundaries.  The  example  provided  based  on  an  international 
collaboration  between  three  towing  tank  facilities  for  purposes  of  procuring  benchmark 
CFD  validation  data  for  ship  hydrodynamics  resistance  and  propulsion  geometry  and 
conditions.  Overlapping  tests  conducted  for  evaluation  of  facilities;  measurement 
systems;  test  procedures;  uncertainty  assessments;  model  size,  offsets,  and  turbulence 
stimulation;  and  facility/model  geometry  and  scale  effect  biases.  The  facilities  were 
David  Taylor  Model  Basin  (DTMB),  Bethesda,  MD  USA;  Istituto  Nazionale  per  Studi  ed 
Esperienze  di  Architettura  Navale  (INSEAN),  Rome,  Italy;  and  Iowa  Institute  of 
Hydraulic  Research  (IIHR),  Iowa  City,  IA,  USA.  Hereafter  designated  as  facilities  A,  B, 
and  C,  respectively.  The  model  geometry  is  DTMB  surface  combatant  5415.  Between  all 
three  facilities,  many  conditions  and  physics  are  under  investigation.  The  data  used  as 
one  of  three  test  cases  at  the  recent  Gothenburg  2000  Workshop  on  Numerical  Ship 
Hydrodynamics  [15].  DTMB  5415  conceived  by  USA  Navy  as  a  preliminary  design  for 
a  surface  combatant  ca.  1980  with  a  sonar  dome  bow  and  transom  stern.  A  and  B  used 
5.72  m,  1/24.8  scale  models,  whereas  C  used  a  3.048  m,  1/46.6  scale  model,  as  shown  in 
Fig.  3.  Thus,  scale  effects  (subgroup  analysis)  considered  in  comparing  C  with  A  and  B. 
The  uncertainty  assessment  procedures  closely  follow  [13]  recommendations  based  on 
[6-7],  Ref.  [16]  provides  an  overview  of  the  overall  results  of  the  collaboration. 
Although  number  of  facilities  is  a  minimum,  the  results  demonstrate  usefulness  of 
approach.  Since  N<10,  Px  is  estimated  using  equation  (13)  as  opposed  to  (12). 
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3.1.  Overlapping  Test  Design,  Data  Reduction  Equations,  Conditions,  Procedures, 
and  Uncertainty  Analysis 

The  most  typical  towing-tank  tests  selected  for  the  overlapping  tests,  i.e.,  resistance, 
sinkage  and  trim,  wave  profile,  wave  elevations,  and  nominal  wake.  Each  institute 
followed  their  usual  procedures;  however,  special  consideration  given  to  integration  of 
uncertainty  assessment  into  all  phases  of  the  experimental  process,  CFD  validation,  and 
complementary  CFD.  Data-reduction  equations  defined  for  residuary  resistance  Cr, 
sinkage  a  and  trim  x,  wave  profile  and  elevations  and  nominal  wake  mean  velocity  V 
and  pressure  Cp.  Similar  conditions  and  locations  were  used  at  all  three  facilities:  Froude 
number  (Fr)  ranges;  Fr  and  spatial  locations  for  uncertainty  analysis;  and  spatial 
resolution  for  wave  elevations  and  nominal  wake. 

Initial  analysis  of  the  results  and  attempt  at  identification  of  facility  biases  done  by 
[16]  using  comparisons  of  differences  between  facilities  (A-B,  A-C,  B-C)  and  the  RSS  of 
their  uncertainties.  However,  this  approach  lacks  a  reference  value  such  that  the 
estimated  facility  biases  depend  on  which  facilities  compared.  Subsequently  all  results 
were  reanalyzed  according  to  the  present  approach.  Presentation  of  the  results  for  all 
tests  is  extensive  and  not  necessary  for  the  purpose  of  demonstration  of  the  usefulness  of 
the  present  approach.  Herein,  the  results  for  the  resistance  test  presented.  The  results  for 
residuary  resistance  follow,  whereas  the  results  for  other  variables  are  included  in 
Appendix  B  and  C. 

The  data  reduction  equation  for  the  resistance  test  is 


Cr  =  Cjm -Clm(l  +  k) 

(35) 

nTm  g 

L  T  -  i 

0.5pU2cS 

(36) 

r.  0.075 

CF  ~  7 

(37) 

(log  10  Re- 2.0) 


The  residuary  resistance  Cr  used  since  it  approximately  removes  a  portion  of  Reynolds 

number  (Re)  scale  effects  due  to  skin  friction  (but  not  wave  breaking).  Cj"'is  the 

frictional  resistance  at  the  measured  towing  tank  temperature  Tm,  k  is  the  form  factor,  Mx 

is  the  force  in  the  axial  direction  (resistance),  p  is  the  towing-tank  water  density,  Uc  is  the 

carriage  speed,  and  S  is  the  design-offsets  wetted  surface  area  for  the  static  condition.  Mx 
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in  kg  is  converted  to  Newtons  by  multiplication  with  g  (gA=9.8009;  gB=9.8033; 
gc=9.8031  m/s')  based  on  the  local  latitude.  Cf  and  k  are  calculated  as  recommended  by 
[17]  using  the  model-ship  correlation  line  equation  (37)  and  Prohaska’s  method, 
respectively. 

Tests  at  A  were  performed  in  basin  no.  2  (575  m  long,  15.5  m  wide,  6.7  m  deep), 
which  is  equipped  with  an  electro-hydraulically  operated  drive  carriage  and  capable  of 
speeds  of  10.3  m/s.  Sidewall  and  end  wall  beaches  enable  20-minute  intervals  between 
carriage  runs.  Washington  Suburban  Sanitation  Commission  supplied  the  towing-tank 
water.  Tests  at  B  were  performed  in  towing  tank  no.  2  (220  m  long,  9  m  wide,  3.6  m 
deep),  which  is  equipped  with  a  single  drive  carriage  that  is  capable  of  speeds  of  10  m/s. 
Sidewall  and  end  wall  beaches  enable  20-minute  intervals  between  carriage  runs.  Natural 
springs  supplied  the  towing  tank  water.  Tests  at  C  were  performed  in  the  IIHR  towing 
tank  (100  m  long  and  3.048  m  wide  and  deep),  which  is  equipped  with  an  electric-motor 
operated  drive  carriage  that  is  cable  driven  by  a  15-horsepower  motor  and  capable  of 
speeds  of  3  m/s.  Sidewall  and  end  wall  beaches  enable  twelve-minute  intervals  between 
carriage  runs.  City  of  Iowa  City  supplied  the  towing  tank  water. 

Equation  (36)  consists  of  individual  MS  for  resistance,  density,  carriage  speed,  and 
surface  area.  Resistance  measured  using  load  cells  and  PC  data  acquisition  and 
reduction,  including  statistical  analysis  of  the  sample  population  (average,  standard 
deviation,  minimums,  maximums,  outliers).  Outliers  identified  and  deleted  using 
Chauvenet’s  criterion.  Density  is  determined  from  measured  Tm  using  fresh  water 
values  as  recommended  by  ITTC  Quality  Manual  Procedure  4.9-03-01-03  Density  and 
Viscosity  of  Water.  Tm  is  measured  daily  using  thermometers.  Carriage  speed  measured 
using  encoders  and  PC  data  acquisition  and  reduction.  Surface  area  is  measured  using 
templates  for  estimating  accuracy  of  model  offsets  and  rulers  and  weights  for  estimating 
accuracy  of  installation  depending  on  whether  model  installed  by  waterline  or 
displacement. 

A  used  a  variable  reluctance,  in-house  manufactured  load  cell,  signal  conditioner,  and 
16-bit  AD  card.  The  load  cell,  signal  conditioner,  and  carriage  PC  AD  card  statically 
calibrated  on  a  test  stand  to  determine  the  voltage-mass  relationship.  Data  acquired  by 
collection  of  2000  discrete  samples  over  5  seconds  at  400  Hz.  Data  filtered  through  a  10 
Hz  low-pass  filter.  B  used  a  Hottinger  Baldwin  Messtechnik  model  Ul,  50  kg  load  cell, 
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signal  conditioner,  and  16-bit  AD  card.  The  load  cell,  signal  conditioner,  and  AD  card 
statically  calibrated  on  a  Kempf  and  Remmers  precision  test  stand  to  determine  the 
voltage-mass  relationship.  Data  acquired  by  collection  of  300  discrete  samples  over  10 
seconds  at  30  Hz.  Amplified  analog  voltages  converted  to  frequency  (3000±2500Hz)  for 
transmission  to  the  AD  card  to  reduce  signal  sensitivity  to  noise.  Data  filtered  through  a 
10  Hz  low-pass  filter.  C  used  a  Nisshio  strain-gage  type  20  kg  load  cell,  signal 
conditioner,  and  12-bit  AD  card.  The  load  cell,  signal  conditioner,  and  AD  card  statically 
calibrated  on  an  IIHR  test  stand  to  determine  the  voltage -mass  relationship.  Data 
acquired  by  collection  of  2000  discrete  samples  over  10  seconds  at  200  Hz.  Data  filtered 
through  a  3  Hz  low-pass  filter. 

Uncertainty  in  Cr  is  equivalent  to  uncertainty  i n  C]"' ;  since,  uncertainties  in  k  and  Cf 
not  considered.  Bias  limits  were  estimated  for  individual  resistance,  density,  carriage 
speed,  and  surface  area  MS,  whereas  precision  limits  were  estimated  end-to-end  using 
equation  (36).  Table  1  summarizes  the  calibration,  data  acquisition,  and  data-reduction 
bias  limits  considered  for  each  MS.  Precision  limits  conducted  over  a  time-period  during 
which  test  conditions  varied  and  in  some  cases  including  reinstallation  of  the  model. 

4.  RESULTS 

4.1.  N-order  level 

Fig.  4  compares  the  individual  facility  (A,  B,  C)  results  for  range  of  Fr.  All  three 
facilities  conducted  uncertainty  analysis  for  Fr=  0.1,  0.28,  and  0.41.  Table  2  provides  N- 
order  level  residuary  resistance  values,  bias  and  precision  limits,  and  total  uncertainties. 
Trends  for  all  three  facilities  are  similar.  Bias  limits  predominate  for  all  Fr.  Although  not 
included  in  Table  2,  bias  limits  for  resistance  and  especially  carriage  speed  are  large  for 
all  speeds,  whereas  bias  limit  for  surface  area  only  significant  for  large  speed.  Precision 
limits  increase  for  increasing  Fr.  Total  uncertainty  decreases  for  increasing  Fr. 

4.2  MxN -order  level 

In  addition,  Table  2  provides  MxN-order  level  facility  uncertainties  and  mean- 
facility  bias  and  precision  limits  and  total  uncertainties  and  facility  certification  or  biases. 
Even  for  M=3,  U x  >  U-  so  UD  ~  U x  .  For  low  speed,  D.  «  UD  and  all  three  facilities 
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are  certified  albeit  at  a  large  interval  (average  17.9%).  Reduction  interval  of  certification 
largely  requires  reduction  individual  facility  biases  for  resistance  and  carriage  speed.  For 
medium  speed,  Dj  >UD  and  all  three  facilities  have  facility  biases  (average  2.3%)  and 

total  uncertainties  larger  than  individual  facility  estimates,  especially  for  C.  For  high 
speed,  D.  >  UD  and  facility  biases  and  total  uncertainties  for  B  and  C  are  large 

(uncertainty  analysis  not  available  for  A).  Fig  4  and  Table  2  suggest  scale  effects 
important  for  C,  for  medium  and  high  speed,  as  shown  next  using  subgroup  analysis. 

4. 3  Subgroup  Analysis 

Table  3  is  similar  to  Table  2,  but  for  subgroup  analysis  in  order  to  isolate  scale 
effects  due  to  smaller  model  size  used  at  C.  In  this  case  mean  facility  based  only  on 
facilities  A  and  B,  which  used  same  model  size,  and  subgroup  analysis  based  on  facility 
C.  For  low  speed,  scale  effects  not  discemable  and  facility  certification  similar  as  before. 
For  medium  and  high  speed,  conclusions  different  from  before,  i.e.,  for  medium  speed 
facilities  A  and  B  certified  at  about  2%  interval,  whereas  facility  C  has  nearly  8%  interval 
facility  bias.  For  high  speed,  facility  B  has  small  and  facility  C  large  facility  biases. 

5  SUMMARY  AND  CONCLUSIONS 

A  statistical  approach  is  set  forth  for  assessing  probabilistic  confidence  intervals 
(i.e.,  intervals  of  certification  or  biases  of  facilities  of  MS)  for  facilities  or  MS  for  specific 
tests,  data  reduction  equations,  conditions,  procedures,  and  uncertainty  analysis  based  on 
MxN-order  level  testing  and  use  of  mean  facility  or  MS  as  reference  values.  MxN-order 
level  testing  defined  as  M  repetitions  of  the  same  N-order  level  experiment  in  M  different 
facilities  or  in  the  same  facility  with  M  different  measurement  systems.  Similarly, 
subgroup  analysis  performed  for  isolating  and  assessing  levels  of  differences  due  to  use 
of  different  model  sizes  (scale  effects)  or  measurement  systems. 

An  example  provided  for  towing  tank  facilities  for  resistance  tests  using  standard 
uncertainty  analysis  procedures  based  on  an  international  collaboration  between  three 
facilities:  two  using  larger  models  and  one  using  smaller  model.  Although  the  number  of 
facilities  is  a  minimum,  the  results  demonstrate  usefulness  of  approach,  including 
subgroup  analysis  for  isolating  differences  due  to  use  of  different  model  sizes.  For  low 
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speed,  all  three  facilities  are  certified,  but  at  a  large  interval  (average  17.9%).  Reduction 
interval  of  certification  largely  requires  reduction  individual  facility  biases  for  resistance 
and  carriage  speed.  For  medium  speed,  facilities  with  larger  models  certified  at  about 
2%,  whereas  facility  with  smaller  model  shows  7.2%  facility  bias.  For  high  speed, 
facilities  with  larger  and  smaller  models  show  2.9%  and  9.3%  facility  biases, 
respectively.  For  certified  facilities,  interval  certification  provides  additional  confidence 
in  accuracy  measurements.  For  non-certified  facilities,  accounting  for  facility  biases 
provides  improved  individual  facility  uncertainties.  Presumably,  design  sets  the 
requirements  on  appropriate  intervals  for  certification  or  biases  of  facilities.  The  results 
are  also  useful  for  developing  strategies  for  reduction  intervals  of  individual  facility 
uncertainties  and  certification  or  facility  biases. 

It  is  reasonable  to  expect  that  results  from  more  facilities  for  same  model 
geometry  with  uncertainty  analysis  will  in  fact  provide  confirmation  of  normal 
distributions  for  MxN-order  level  testing  and  improved  estimates  for  intervals  of 
certification  or  biases  of  facilities.  This  based  on  previous  work  between  more  facilities 
using  same  model  geometry  but  without  uncertainty  analysis  and  using  uncertainty 
analysis  but  with  different  model  geometry.  It  is  also  reasonable  to  expect  that  situation 
will  be  similar  for  different  types  of  facilities  and  MS.  Both  these  expectations  support 
recommendation  of  future  collaborations  between  more  facilities.  International 
collaborations  are  attractive  from  resource  perspective  and  in  achieving  ground  truth. 
Knowledge  of  intervals  of  certification  or  biases  is  important  for  design,  accrediting 
facilities  or  measurement  systems,  and  CFD  validation. 
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TABLE  AND  FIGURES 


Table  1  Bias  limits  for  resistance,  density,  carriage  speed,  and  surface  area 


Calibration 

Data  Acquisition 

Data  Reduction 

Resistance 

Load  Cell 

Curve  fit 

Alignment  load  cell 

Density 

Thermometer 

Carriage  Speed 

Encoder 

Wheel  diameter 

Time  base 

AD  conversions 

Curve  fit 

Surface  area 

Template 

Rulers 

Weights 

Table  2.  N-order  level  residuary  resistance  values,  bias  and  precision  limits,  and  total 
uncertainties;  MxN-order  level  facility  uncertainties  and  mean-facility  bias  and  precision 
limits  and  total  uncertainties;  and  facility  certification  or  biases. 


Fr  & 


Facility  Certification 

N-order  level  (%  Xi ) 

MxN-order  level  (%  X  ) 

or 

biases  (%  X  ) 

Xt  B2Xi  Pi  UXi 

Ux< 

4  4  ux 

A 

UD, 

uFBi 

uTi 

6.00E-04  76.3%  23.7%  10.4% 

11.3% 

8.6% 

14.5% 

0 

11.3% 

5.23E-04  69.4%  30.6%  21.2% 

20.1% 

75.7%  24.3%  9.1% 

-5.3% 

22.0% 

0 

20.1% 

5.34E-04  87.6%  12.4%  14.9% 

14.5% 

-3.3% 

17.1% 

0 

14.5% 

5.52E-04  77.8%  22.2%  15.5% 

15.3% 

0.0% 

17.9% 

0 

15.3% 

1.33E-03  45.5%  54.5%  1.1% 

1.1% 

3.0% 

1.6% 

2.6% 

2.8% 

1.32E-03  80.0%  20.0%  2.1% 

2.1% 

81.6%  18.4%  1.2% 

2.2% 

2.4% 

0 

2.1% 

1.22E-03  89.2%  10.8%  2.7% 

2.6% 

-5.3% 

2.8% 

4.4% 

5.1% 

1.29E-03  71.6%  28.4%  2.0% 

1.9% 

0.0% 

2.3% 

2.3% 

3.4% 

3.69E-03  NA  NA  NA 

NA 

-0.2% 

NA 

NA 

NA 

3.94E-03  66.3%  33.7%  1.1% 

1.2% 

73.6%  26.4%  0.8% 

6.7% 

1.4% 

6.5% 

6.6% 

3.45E-03  80.5%  19.5%  1.3% 

1.2% 

-6.5% 

1.4% 

6.3% 

6.5% 

3.69E-03  73.4%  26.6%  1.2% 

1.2% 

0.0% 

1.4% 

6.4% 

6.5% 
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Table  3.  Subgroup  Analysis:  N-order  level  residuary  resistance  values,  bias  and  precision 
limits,  and  total  uncertainties;  MxN-order  level  facility  uncertainties  and  mean  facility 
bias  and  precision  limits  and  total  uncertainties;  and  facility  certification  or  biases 


Fr  & 

Facility 

N-order  level  (%  Xt ) 

MxN-order  level  (%  X  ) 

Facility  Certification 

or  biases  (%  X  ) 

X,  Pi  UXj 

uXi 

B-  P-  U- 

X  X  ^  X 

A 

uDl 

uFBi 

uTi 

0.10 

A 

6.00E-04  76.3%  23.7%  10.4% 

11.1% 

71.1%  28.9%  11.3% 

6.8% 

15.9% 

0 

11.1% 

B 

5.23E-04  69.4%  30.6%  21.2% 

19.8% 

-6.8% 

22.8% 

0 

19.8% 

C 

5.34E-04  87.6%  12.4%  14.9% 

14.2% 

-4.8% 

18.2% 

0 

14.2% 

AVE 

5.61E-04  77.8%  22.2%  15.5% 

15.0% 

-1.6% 

18.9% 

0 

15.0% 

0.28 

A 

1.33E-03  45.5%  54.5%  1.1% 

1.1% 

72.8%  27.2%  1.2% 

0.4% 

1.6% 

0 

1.1% 

B 

1.32E-03  80.0%  20.0%  2.1% 

2.1% 

-0.4% 

2.4% 

0 

2.1% 

C 

1.22E-03  89.2%  10.8%  2.7% 

2.5% 

-7.7% 

2.8% 

7.2% 

7.6% 

AVE 

1.33E-03  71.6%  28.4%  2.0% 

1.9% 

-2.6% 

2.3% 

2.4% 

3.6% 

0.41 

A 

3.69E-03  NA  NA  NA 

NA 

66.3%  33.7%  1.1% 

-3.3% 

NA 

NA 

NA 

B 

3.94E-03  66.3%  33.7%  1.1% 

1.1% 

3.3% 

1.6% 

2.9% 

3.1% 

C 

3.45E-03  80.5%  19.5%  1.3% 

1.2% 

-9.5% 

1.6% 

9.3% 

9.4% 

AVE 

3.81E-03  73.4%  26.6%  1.2% 

1.1% 

-3.2% 

1.6% 

6.1% 

6.3% 

Figure  1.  MxN-order  level  testing 
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Fig  2  Estimating  intervals  of  certification  or  biases  of  facilities  using  mean  as  reference  value: 
facility  i  certified  at  UD  and  facility  j  with  facility  bias  UFB  . 


(a)DTMB  model  5415 


(b)  INSEAN  model  2340A 


(c)  IIHR  model  5512 
Fig  3  Surface  combatant  model  ships 
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0.006 


Figure  4  Individual  facility  residuary  resistance  results  and  mean  facility  residuary 
resistance  and  individual  facility  uncertainty  bands  at  Fr=0.1,  0.28  and  0.41. 
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APPENDIX  A:  ESTIMATING  INTERVALS  OF  CERTIFICATION  OR  BIASES 
OF  FACILITIES  OR  MS  USING  REFERENCE  VALUES 

Assume  a  reference  value  known  for  the  experimental  result  designated  XR  with 
uncertainty  UR  both  of  which  considered  as  standard  values.  Although  UR  is  likely  much 
less  than  U x  or  U-,  it  is  retained  for  completeness. 

For  the  mean  facility,!)  is  defined  as  the  difference  between  the  mean  facility X 
and  reference  XR  values 

D=X-Xr  (Al) 

and  its  uncertainty  U-  is  defined  as  the  RSS  of  the  uncertainties  of  X  and  XR 

U-D  =  V  uf+ul  (A2) 

If  the  absolute  value  of  D  is  less  than  U- 

\D\<Un  (A3) 

then  the  mean  facility  is  certified  at  the  interval  O'  - ,  whereas  if  the  absolute  value  of  D 
is  greater  than  U  r) 

\D\>UB  (A5) 


then  the  mean  facility  bias  U  ^  is  defined  as 


U—  =D2 -U- 

^  FB  ^  D 


with  total  uncertainty  U T 


Uf=U2,+U 


2 

FB 


(A6) 


(A7) 


If  D  is  much  greater  than  U  f) 


D  »UD 


(A8) 
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then  D  approximately  equals  the  mean  facility  bias  error  /3FB 

D  =  PFB  (A9) 

such  that  the  biases  can  be  estimated  in  both  sign  and  magnitude  and  used  for  calibration. 
For  the  individual  facility,!),  is  defined  as  the  difference  between  the  individual 

facility  Xt  and  reference  XR  values 

Di=Xi-XR  (A10) 

and  its  uncertainty  UD  is  defined  as  the  RSS  of  the  uncertainties  of  A)  and  XR 

un=j  Ul+Ul  (All) 

If  the  absolute  value  of  Z)  is  less  than  UD 

|  Di\<UDi  (A  12) 

then  the  individual  facility  is  certified  at  the  interval  UD  ,  whereas  if  the  absolute  value  of 
Z)  is  greater  than  UD 

\Di\>UDi  (A  13) 

then  the  individual  facility  bias  UFB  is  defined  as 

U2FBi  =  Df -U2Dt  (A  14) 

with  total  uncertainty  UTj 

U2T=U2x+U2FBi  (A  15) 

If  |Z);  |  is  much  greater  than  U D 

\Di\»UDi  (A  16) 

then  Dj  approximately  equals  the  individual  facility  bias  error  fiFB 

D  =  PFB<  (A  17) 


such  that  the  biases  can  be  estimated  in  both  sign  and  magnitude  and  used  for  calibration. 
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APPENDIX  B:  FACILITY  CERTIFICATION/BIASES  FOR  SINKAGE,  TRIM, 
WAVE  PROFILE,  WAVE  ELEVATION,  AND  NOMINAL  WAKE  VELOCITY 

USING  AVERAGE  X  =  (A+B+C)/3 


Note  for  Table  2  in  text: 

(1)  Xf  for  A  and  C  from  Interpolated  data  files  (dtmbintp.tec,  iihr  intp.tec) 

(2)  A  for  B  from  IIHR  Report  421  (pi 8)  for  INSEAN 

(3)  All  B\  and  P]  from  [16] 

(4)  U x  for  A  from  DTMB’s  Raw  data  files  (dtmb_resist.tec) 

(5)  U x  for  B  from  Angelo  Olivieri’s  email  to  Fred  Stern  (04/08/2003) 

(6)  U x  for  C  from  IIHR  Raw  data  document 
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Table  Bl.  Facility  certification  or  biases  for  sinkage  results  cr 


Fr  & 

Facility 

N-order  level  (%  A.) 

MxN-order  level  (%  X  ) 

Facility  Certification 

or  biases  (%  X  ) 

Xt  Bl  Pi  Ux_ 

uXl 

4  4  ux 

D, 

uDi 

uFBi 

uTi 

0.10 

A 

7.49E-04  75.6%  24.4%  12.2% 

13.4% 

8.8%  91.2%  15.3% 

10.1% 

20.3% 

0 

13.4% 

B 

6.98E-04  0.0%  100.0%  42.0% 

43.1% 

2.6% 

45.7% 

0 

43.1% 

C 

5.93E-04  82.2%  17.8%  8.7% 

7.6% 

-12.8% 

17.1% 

0 

7.6% 

AVE 

6.80E-04  52.6%  47.4%  21.0% 

21.4% 

0.0% 

27.7% 

0 

21.4% 

0.28 

A 

7.35E-03  68.4%  32.6%  5.6% 

5.5% 

39.3%  60.7%  2.5% 

-1.0% 

6.1% 

0 

5.5% 

B 

7.39E-03  0.0%  100.0%  4.7% 

4.7% 

-0.3% 

5.3% 

0 

4.7% 

C 

7.51E-03  30.4%  69.6%  1.4% 

1.4% 

1.3% 

2.9% 

0 

1.4% 

AVE 

7.42E-03  32.9%  67.4%  3.9% 

3.9% 

0.0% 

4.8% 

0 

3.9% 

0.41 

A 

1.73E-02  56.3%  44.7%  2.5% 

2.4% 

21.9%  78.1%  1.3% 

-5.1% 

2.7% 

4.3% 

4.9% 

B 

1.88E-02  0.0%  100.0%  2.9% 

3.0% 

3.3% 

3.3% 

0 

3.0% 

C 

1.85E-02  42.8%  57.2%  0.6% 

0.6% 

1.8% 

1.4% 

1.1% 

1.2% 

AVE 

1.82E-02  33.0%  67.3%  2.0% 

2.0% 

0.0% 

2.5% 

1.8% 

3.1% 

Note: 

(1)  DRE  for  sinkage  in  ONR  paper  figures  for  three  institutes  is 

2(AFP  +  AAP)  .  ,  „  2  2(AFP  +  AAP) 

cr  = -  instead  ol  cr  =  — - - 

L  Fr 2  2  L 

(2)  Uncertainties  in  terms  of  percentages  remain  the  same  for  both  above  DREs,  while 
uncertainty  magnitudes  are  different  as  long  as  we  assume  that  Fr  number  has  no 
uncertainty  for  simplicity. 

(3)  Xi  for  A,  B,  and  C  from  Interpolated  data  files  (dtmbintp.tec,  sean  intp.tec, 
iihrintp.tec) 

(4)  All  B\  and  Pi  from  [16] 

(5) B\  and  Px  for  DTMB  (Fr=0.28  and  0.41)  are  best  estimates  from  [16] 

(6)  All  U x  based  on  mean  values  from  Interpolated  data  files  (dtmb  intp.tec, 
sean  intp.tec,  iihr  intp.tec) 
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Table  B2.  Facility  certification  or  biases  for  trim  results  z 


Fr  & 

Facility 

N-order  level  (%  Xt) 

MxN-order  level  (%  X  ) 

Facility  Certification 

or  biases  (%  X  ) 

X ,  B%  Pi  UXi 

A, 

4  p}  u-x 

A 

A, 

UFB, 

uTj 

0.10 

A 

-2.05E-04  64.5%  35.5%  14.4% 

6.1% 

6.0%  94.0%  14.5% 

-57.4% 

15.7% 

55.2% 

55.5% 

B 

-6.13E-04  0.0%  100.0%  32.0% 

40.9% 

27.7% 

43.3% 

0 

40.9% 

C 

-6.22E-04  50.8%  49.2%  10.2% 

13.3% 

29.7% 

19.6% 

22.3% 

25.9% 

AVE 

-4.80E-04  38.4%  61.6%  18.9% 

20.1% 

0.0% 

26.2% 

25.8% 

40.8% 

0.28 

A 

-3.90E-03  54.7%  46.3%  2.8% 

2.6% 

18.1%  81.9%  1.8% 

-5.8% 

3.2% 

4.8% 

5.5% 

B 

-3.77E-03  0.0%  100.0%  4.7% 

4.3% 

-8.9% 

4.7% 

7.6% 

8.8% 

C 

-4.75E-03  36.1%  63.9%  1.8% 

2.1% 

14.7% 

2.8% 

14.5% 

14.6% 

AVE 

-4.14E-03  30.3%  70.1%  3.1% 

3.0% 

0.0% 

3.5% 

9.0% 

9.6% 

0.41 

A 

1.36E-02  38.1%  61.9%  1.5% 

1.6% 

18.9%  81.1%  0.8% 

6.6% 

1.8% 

6.4% 

6.6% 

B 

1.40E-02  0.0%  100.0%  0.9% 

1.0% 

10.1% 

1.2% 

10.0% 

10.0% 

C 

1.06E-02  4.1%  95.9%  1.8% 

1.5% 

-16.7% 

1.7% 

16.6% 

16.6% 

AVE 

1.27E-02  14.1%  85.9%  1.4% 

1.3% 

0.0% 

1.6% 

11.0% 

11.1% 

Note: 

(1)  DRE  for  trim  in  ONR  paper  figures  for  three  institutes  is 

r = instead of r = 

L  Fr 2  L 

(2)  Uncertainties  in  terms  of  percentages  remain  the  same  for  both  above  DREs,  while 
uncertainty  magnitudes  are  different  as  long  as  we  assume  that  Fr  number  has  no 
uncertainty  for  simplicity. 

(3)  Xt  for  A,  B,  and  C  from  Interpolated  data  files  (dtmbintp.tec,  sean  intp.tec, 
iihrintp.tec) 

(4)  All  B\  and  P]  from  [16] 

(5 ) BX  and  Py  for  DTMB  (Fr=0.28)  are  best  estimates  from  [16] 

(6)  All  U x  based  on  mean  values  from  Interpolated  data  files  (dtmb  intp.tec, 
sean  intp.tec,  iihr  intp.tec) 
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Table  B3.  Facility  certification  or  biases  for  wave  profile  results  C, 


Fr& 

Facility 

N-order  level  (%  Dx  ) 

MxN-order  level  (%  D  v  ) 

Facility  Certification 

or  biases  (%  D  x  ) 

X,  DXi  B\  P;  UXi 

uXi 

B 2  P-  U- 

X  X  ^  X 

A 

uDl 

UFB, 

UTi 

0.28 

A 

1.82E-02  1.89E-02  64.5%  35.5%  3.5% 

3.2% 

85.9%  14.1%  2.2% 

-6.4% 

3.9% 

5.1% 

6.0% 

B 

2.00E-02  2.09E-02  100.0%  0.0%  4.2% 

4.2% 

2.4% 

4.7% 

0 

4.2% 

C 

2.03E-02  2.22E-02  83.7%  16.3%  3.4% 

3.7% 

3.9% 

4.3% 

0 

3.7% 

AVE 

1.95E-02  2.07E-02  82.7%  17.3%  3.7% 

3.7% 

0.0% 

4.3% 

1.7% 

4.6% 

0.41 

A 

3.02E-02  3.62E-02  64.5%  35.5%  1.8% 

2.0% 

81.5%  18.5%  1.2% 

2.0% 

2.3% 

0 

2.0% 

B 

2.82E-02  2.46E-02  100.0%  0.0%  2.6% 

1.9% 

-4.2% 

2.3% 

3.5% 

4.0% 

C 

3.03E-02  3.86E-02  81.6%  18.4%  2.0% 

2.3% 

2.2% 

2.6% 

0 

2.3% 

AVE 

2.95E-02  3.32E-02  82.0%  18.0%  2.1% 

2.1% 

0.0% 

2.4% 

1.2% 

2.8% 

A  is  the  maximum  elevation  on  the  wave  profile. 


Note: 

(1)  Xi  from  Data  files  before  interpolation  (dtmb_dpwpro28.tec,  sean_dpwpro28.tec, 
iihr_dpwpro28.tec) 

(2)  Dx  for  A  from  DTMB’s  website:  http://www50.dt.navv.mil/54 1 5/profilc.hlml,  (0.15 
inch  for  all  points) 

(3)  Dx  for  B  from  IIHR  Report  421  (p.53)  for  INSEAN 

(4)  Dx  for  C  from  IIHR’s  raw  data  files  (wpro2801.tec,  wpro2802.tec,  wpro2803.tec, 
wpro4101.tec,  wpro4102.tec,  wpro4103.tec) 

(5)  All  BXj  and  PXj  from  [16] 

(6)  All  UXj  from  [16] 
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Table  B4.  Facility  certification  or  biases  for  wave  elevation  results  C,  at  cut  y=0.324 


Fr  & 

Facility 

N-order  level  (%  Dx  ) 

MxN-order  level  (%  Dx  ) 

Facility  Certification 

or  biases  (%  Dx  ) 

X,.  DXt  B\  Pi  UXt 

A, 

4  4  ux 

A 

uDt 

UFB, 

UTj 

0.28 

4.19E-03  1.17E-02  76.6%  23.4%  2.7% 

2.8% 

66.0%  34.0%  1.7% 

-1.8% 

3.3% 

0 

2.8% 

B 

5.20E-03  1.09E-02  64.9%  35.1%  2.4% 

2.3% 

7.2% 

2.9% 

6.6% 

7.0% 

3.79E-03  1.11E-02  59.0%  41.0%  3.4% 

3.4% 

-5.4% 

3.8% 

3.9% 

5.1% 

4.39E-03  1.12E-02  66.8%  33.2%  2.9% 

2.9% 

0.0% 

3.3% 

3.5% 

5.0% 

Xt  is  the  maximum  elevation  on  the  cut  y=0.324. 


Note: 

(1)  X.  from  Interpolated  data  files  (dtmb_int324.tec,  insean_int324.tec,  iihr_int324.tec) 

(2)  Dx  for  A  from  Interpolated  data  files  (dtmbintpat.tec) 

(3)  Dx  for  B  from  IIHR  Report  421  (p.54)  for  INSEAN 

(4)  Dx  for  C  from  (Joe  Longo)  IIHR’s  raw  data  file  (a0_101.dat  zone  =  “Steady”) 

(5)  All  B2x  and  Px  from  [16] 

(6)  All  UXj  from  [16] 

(7 ) BX  and  Px  for  DTMB  (Fr=0.28)  are  best  estimates  from  [16] 


27 


Table  B5.  Facility  certification  or  biases  for  wave  elevation  results  C, 


Fr  & 

Facility 

N-order  level  (%  Dx  ) 

MxN-order  level  (%DV  ) 

Facility  Certification 

or  biases  (%  Dx  ) 

DXl 

P2X,  Ux> 

uXl 

4  p=  ux 

Di 

UD, 

uFBt 

UTj 

0.28 

A 

1.17E-02 

76.6%  23.4%  2.7% 

2.8% 

66.0%  34.0%  1.7% 

4.1%  3.3%  2.5%  3.8% 

B 

1.09E-02 

64.9%  35.1%  2.4% 

2.3% 

2.8%  2.9%  0  2.3% 

C 

1.11E-02 

59.0%  41.0%  3.4% 

3.4% 

2.0%  3.8%  0  3.4% 

AVE 

1.12E-02 

66.8%  33.2%  2.9% 

2.9% 

2.9%  3.3%  0.8%  3.2% 

Note: 

(1)  Dx  for  A  from  Interpolated  data  files  (dtmbintpat.tec) 

(2)  Dx  for  B  from  IIHR  Report  421  (p.54)  for  INSEAN 

(3)  Dx  for  C  from  (Joe  Longo)  IIHR’s  raw  data  fdes  (a0_101.dat  zone  =  “Steady”) 

(4)  All  B\  and  Px  from  [16] 

(5)  All  UXj  from  [16] 

(6 ) BX  and  Px  for  DTMB  (Fr=0.28)  are  best  estimates  from  [16] 

(7) A11  D.  calculated  from  Interpolated  data  files  (dtmb  intpat.tec,  insean  intpat.tec, 
iihrintpat.tec) 
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Table  B6.  Facility  certification  or  biases  for  nominal  wake  velocity  results  at  Fr=0.28 


N-order  level  (%  Dx  ) 

Facility  Certification 

Fr  & 

MxN-order  level  (%  Dx  ) 

or 

biases  (%  D 

xt) 

Facility 

LTR 

HTR 

SQRT  (LTRA2+HTRA2)/2 

A, 

B 2  P2 

Xi  Xi 

Ux, 

B 2  P2 

Xi  Xi 

uXl 

Bl 

P 2 

Xi 

uXl 

Ux, 

B-  P- 

X  1  X 

Ux 

A 

uDi 

uFBi 

uTt 

A 

1.00E+00 

65.4%  34.4% 

1.6% 

74.5%  26.5% 

12.5% 

73.6% 

26.4% 

6.3% 

6.3% 

1.5% 

6.6% 

0 

6.3% 

0.28 

B 

1.02E+00 

47.8%  52.2% 

0.4% 

60.4%  39.6% 

0.4% 

53.5% 

46.5% 

0.3% 

0.3% 

75.3%  24.7% 

2.2% 

1.3% 

2.2% 

0 

0.3% 

U 

C 

1.01E+00 

99.8%  0.2% 

1.2% 

99.2%  0.8% 

3.1% 

99.3% 

0.7% 

1.7% 

1.7% 

1.7% 

2.7% 

0 

1.7% 

AVE 

1.01E+00 

71.0%  28.9% 

1.1% 

78.0%  22.3% 

5.3% 

75.5% 

24.5% 

2.8% 

2.7% 

1.5% 

3.9% 

0 

2.7% 

A 

9.62E-02 

54.3%  46.7% 

2.9% 

43.5%  56.5% 

6.5% 

45.2% 

54.8% 

3.6% 

3.0% 

4.1% 

3.4% 

2.2% 

3.7% 

0.28 

B 

1.19E-01 

21.2%  78.8% 

1.9% 

15.9%  84.1% 

2.7% 

17.6% 

82.4% 

1.6% 

1.7% 

71.6%  28.4% 

1.7% 

4.3% 

2.4% 

3.6% 

4.0% 

V 

C 

1.33E-01 

99.7%  0.3% 

5.5% 

93.3%  6.7% 

3.8% 

97.7% 

2.3% 

3.3% 

3.8% 

4.0% 

4.2% 

0 

3.8% 

AVE 

1.16E-01 

58.4%  41.9% 

3.4% 

50.9%  49.1% 

4.3% 

53.5% 

46.5% 

2.9% 

2.8% 

4.1% 

3.3% 

1.9% 

3.8% 

A 

1.47E-01 

65.3%  34.7% 

6.5% 

44.6%  35.4% 

3.7% 

63.3% 

36.7% 

3.6% 

3.3% 

6.0% 

3.6% 

4.8% 

5.8% 

0.28 

B 

1.65E-01 

79.1%  20.9% 

1.0% 

87.9%  12.1% 

0.9% 

83.1% 

16.9% 

0.6% 

0.6% 

81.7%  18.3% 

1.6% 

11.3% 

1.7% 

11.2% 

11.2% 

W 

C 

1.79E-01 

99.9%  0.1% 

4.1% 

99.2%  0.8% 

4.5% 

99.5% 

0.5% 

3.0% 

3.3% 

5.8% 

3.7% 

4.5% 

5.6% 

AVE 

1.64E-01 

81.4%  18.6% 

3.8% 

77.2%  16.1% 

3.0% 

82.0% 

18.0% 

2.4% 

2.4% 

7.7% 

3.0% 

6.8% 

7.5% 

Note: 

(1)  Dx  for  A  from  DTMB  Raw  data  files:  (dtmbdpnwake.tec) 

(2)  Dx  for  B  from  IIHR  Report  421  (p.56)  for  INSEAN,  but  dynamic  range  for  U  is  divided  by  2. 

(3)  Dx  for  C  from  IIHR  Raw  data  document 

(4)  All  B\  and  Px  from  [16] 
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(5)  All  UXj  from  [16] 

(6)  B\  and  Px  for  DTMB  are  best  estimates  from  [16] 

(7)  All  D.  calculated  from  Interpolated  data  files  (dtmb  dpINTnwake.tec,  sean  dpINTnwake.tec,  iihr  dpINTnwake.tec) 
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Figure  B1  Residuary  resistance  results  and  uncertainty  bands  at  Fr=0.1,  0.28  and 

0.41 


Note: 

(1)  DTMB  data  file:  \23rdONRdata\Resistance\dtmb_dpcr.tec 

(2)  INSEAN  data  file:  \23rdONRdata\Resistance\sean_dpcr.tec 

(3)  IIHR  data  file:  \23rdONRdata\Resistance\iihr_dpcr.tec 

(4)  Mean  data  file:  \23rdONRdata\Resistance\meanABC.txt 

(5)  DTMB  uncertainty  file:  \23rdONRdata\Resistance\dtmb_dpcr_un.dat 

(6)  INSEAN  uncertainty  file:  \23rdONRdata\Resistance\sean_dpcr_un.dat 

(7)  IIHR  uncertainty  file:  \23rdONRdata\Resistance\iihr_dpcr_un.dat 
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Figure  B2  Sinkage  results  and  uncertainty  bands  at  Fr=0.1,  0.28  and  0.41 


Note: 

(1)  DTMB  data  file:  \23rdONRdata\Sinkage_Trim\dtmb_dpst.tec 

(2)  INSEAN  data  file:  \23rdONRdata\Sinkage_Trim\sean_dpst.tec 

(3)  IIHR  data  file:  \23rdONRdata\Sinkage_Trim\iihr_dpst.tec 

(4)  Mean  data  file:  \23rdONRdata\Sinkage_Trim\meanABC.txt 

(5)  DTMB  uncertainty  file:  \23rdONRdata\Sinkage_Trim\dtmb_dpcr_sink_un.dat 

(6)  INSEAN  uncertainty  file:  \23rdONRdata\Sinkage_Trim\sean_dpcr_sink_un.dat 

(7)  IIHR  uncertainty  file:  \23rdONRdata\Sinkage_Trim\iihr_dpcr_sink_un.dat 
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Figure  B3  Trim  results  and  uncertainty  bands  at  Fr=0.1,  0.28  and  0.41 


Note: 

(1)  DTMB  data  file:  \23rdONRdata\Sinkage_Trim\dtmb_dpst.tec 

(2)  INSEAN  data  file:  \23rdONRdata\Sinkage_Trim\sean_dpst.tec 

(3)  IIHR  data  file:  \23rdONRdata\Sinkage_Trim\iihr_dpst.tec 

(4)  Mean  data  file:  \23rdONRdata\Sinkage_Trim\meanABC.txt 

(5)  DTMB  uncertainty  file:  \23rdONRdata\Sinkage_Trim\dtmb_dpcr_trim_un.dat 

(6)  INSEAN  uncertainty  file:  \23rdONRdata\Sinkage_Trim\sean_dpcr_trim_un.dat 

(7)  IIHR  uncertainty  file:  \23rdONRdata\Sinkage_Trim\iihr_dpcr_trim_un.dat 
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0.0250 


Figure  B4  Wave  profile  results  and  uncertainty  bands  at  <^max  (Fr=0.28) 


Note: 

(1)  DTMB  data  file:  \23rdONRdata\Wave_Profiles\dtmb_dpwpro28.tec 

(2)  INSEAN  data  file:  \23rdONRdata\Wave_Profiles\sean_dpswpro28.tec 

(3)  IIHR  data  file:  \23rdONRdata\Wave_Profiles\iihr_dpwpro28.tec 

(4)  Mean  data  file:  \23rdONRdata\Wave_Profiles\mean28ABC.txt 

(5)  DTMB  uncertainty  file:  \23rdONRdata\Wave_Profiles\dtmb_dp28_un.dat 

(6)  INSEAN  uncertainty  file:  \23rdONRdata\Wave_Profiles\sean_dp28_un.dat 

(7)  IIHR  uncertainty  file:  \23rdONRdata\Wave_Profiles\iihr_dp28_un.dat 


34 


Figure  B5  Wave  profile  results  and  uncertainty  bands  at  <^max  (Fr=0.41) 


Note: 

(1)  DTMB  data  file:  \23rdONRdata\Wave_Profiles\dtmb_dpwpro41.tec 

(2)  INSEAN  data  file:  \23rdONRdata\Wave_Profiles\sean_dpswpro41.tec 

(3)  IIHR  data  file:  \23rdONRdata\Wave_Profiles\iihr_dpwpro41.tec 

(4)  Mean  data  file:  \23rdONRdata\Wave_Profiles\mean41ABC.txt 

(5)  DTMB  uncertainty  file:  \23rdONRdata\Wave_Profiles\dtmb_dp41_un.dat 

(6)  INSEAN  uncertainty  file:  \23rdONRdata\Wave_Profiles\sean_dp41_un.dat 

(7)  IIHR  uncertainty  file:  \23rdONRdata\Wave_Profiles\iihr_dp41_un.dat 
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Figure  B6  Wave  elevation  results  and  uncertainty  bands  for  cut  y=0.324  (Fr=0.28) 

Note: 

(1)  DTMB  data  file:  \23rdONRdata\Wave_Elevations\dtmb_int324.tec 

(2)  INSEAN  data  file:  \23rdONRdata\Wave_Elevations\sean_int324.tec 

(3)  IIHR  data  file:  \23rdONRdata\Wave_Elevations\iihr_int324.tec 

(4)  Mean  data  file:  \23rdONRdata\Wave_Elevations\meanABC.txt 

(5)  DTMB  uncertainty  file:  \23rdONRdata\Wave_Elevations\dtmb_dp28_un.dat 

(6)  INSEAN  uncertainty  file:  \23rdONRdata\Wave_Elevations\sean_dp28_un.dat 

(7)  IIHR  uncertainty  file:  \23rdONRdata\Wave_Elevations\iihr_dp28_un.dat 
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DTMB  541 5  Far  Field  Wave  Elevation  Results 


-0.05 


INSEAN  2340A  Far  Field  Wave  Elevation  Results 

-0.05  F- 


IIHR  5512  Far  Field  Wave  Elevation  Results 

-0.05  p 


Mean  Far  Field  Wave  Elevation  Results 

-0.05  F- 


Figure  B7  Wave  elevation  results  for  far  field  and  the  mean  result  (Fr=0.28) 


Note: 

(1)  DTMB  data  file:  \23rdONRdata\Wave_Elevations\dtmb_intpat.tec 

(2)  INSEAN  data  file:  \23rdONRdata\Wave_Elevations\insean_intpat.tec 

(3)  IIHR  data  file:  \23rdONRdata\Wave_Elevations\iihr_intpat.tec 

(4)  Mean  data  file:  \23rdONRdata\Wave_Elevations\mean_Zeta_ABC.dat 
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-0.05 
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-0.15 
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-0.35 


D  for  DTMB  5415  Far  Field  Wave  Ele_vation  Results 
Dj  =  4.1%,  UDi  =  3.3%  (%DC) 


a 


D.  for  INSEAN  2340A  Far  Field  Wave  Elevation  Results 
Dj  =  2.8%,  UDi  =  2.9%  (%D;) 


D.  for  IIHR  5512  Far  Field  Wave  Elevation  Results 


Figure  B8  Contours  of  difference  Di  of  wave  elevation  results  (Fr=0.28) 


Note: 

(1)  DTMB  data  file:  \23rdONRdata\Wave_Elevations\DTMB_DZeta_ABC.tec 

(2)  INSEAN  data  file:  \23rdONRdata\Wave_Elevations\SEAN_DZeta_ABC.tec 

(3)  IIHR  data  file:  \23rdONRdata\Wave_Elevations\IIHR_DZeta_ABC.tec 
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Certification  of  Far  Field  Wave  Elevation  Results  (INSEAN) 
Mean  =  (A+B+C)/3 


Certification  of  Far  Field  Wave  Elevation  Results  (IIHR) 
Mean  =  (A+B+C)/3 


Figure  B9  3D  surface  of  difference  Di  of  wave  elevation  results  (Fr=0.28) 

Note: 

(1)  DTMB  data  file:  \23rdONRdata\Wave_Elevations\DTMB_DZeta_ABC.tec 

(2)  INSEAN  data  file:  \23rdONRdata\Wave_Elevations\SEAN_DZeta_ABC.tec 

(3)  IIHR  data  file:  \23rdONRdata\Wave_Elevations\IIHR_DZeta_ABC.tec 
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Figure  BIO  Nominal  wake  velocity  results  and  the  mean  velocity  result  (Fr=0.28) 
Note: 

(1)  DTMB  data  file:  \23rdONRdata\Nominal  Wake\DTMB  dpINTnwake.dat 

(2)  INSEAN  data  fde:  \23rdONRdata\Nominal  Wake\SEAN  dpINTnwake.dat 

(3)  IIHR  data  file:  \23rdONRdata\Nominal  WakeUIHR dpINTnwake.dat 

(4)  Mean  data  file:  \23rdONRdata\Nominal_Wake\MeanABC_UVW.dat 
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Figure  Bll  Contours  of  difference  Di  of  nominal  wake  velocity  results  (Fr=0.28) 
Note: 

(1)  DTMB  data  file:  \23rdONRdata\Nominal_Wake\MeanABC_DA_perc.dat 

(2)  INSEAN  data  fde:  \23rdONRdata\Nominal  Wake\MeanABC_DB  perc.dat 

(3)  IIHR  data  file:  \23rdONRdata\Nominal_Wake\MeanABC_DC_perc.dat 
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Certification  of  nominal  wake  velocity  component  V  (DTM  B) 
Mean  =  (A+B+C)/3 


z 
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Figure  B12  3D  surface  of  difference  Di  of  nominal  wake  velocity  results  (Fr=0.28) 
Note: 

(1)  DTMB  data  fde:  \23rdONRdata\Nominal  Wake\MeanABC_DA.dat 

(2)  INSEAN  data  fde:  \23rdONRdata\Nominal  Wake\MeanABC_DB.dat 

(3)  IIHR  data  file:  \23rdONRdata\Nominal  Wake\MeanABC_DC.dat 
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APPENDIX  C:  FACILITY  CERTIFICATION/BIASES  FOR  SINKAGE,  TRIM, 
WAVE  PROFILE,  WAVE  ELEVATION,  AND  NOMINAL  WAKE  VELOCITY 

USING  AVERAGE  X  =  (A+B)/2 


Note  for  Table  3  in  text: 

(1)  X.  for  A  and  C  from  Interpolated  data  files  (dtmb  intp.tec,  iihr  intp.tec) 

(2)  Xt  for  B  from  IIHR  Report  421  for  INSEAN 

(3)  All  B\  and  P]  from  [16] 

(4)  U x  for  A  from  DTMB’s  Raw  data  files  (dtmb  resist.tec) 

(5)  U x  for  B  from  Angelo  Olivieri’s  email  to  Fred  Stern  (04/08/2003) 

(6)  Ux  for  C  from  IIHR  Raw  data  document 
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Table  Cl.  Facility  certification  or  biases  for  sinkage  results  cr 


Fr  & 

Facility 

N-order  level  (%  Xt) 

MxN-order  level  (%  X  ) 

Facility  Certification 

or  biases  (%  X  ) 

X,  B\  Pi  Ux_ 

uXi 

4  p2x  ux 

D, 

UD, 

uFBi 

UTj 

0.10 

A 

7.49E-04  75.6%  24.4%  12.2% 

12.6% 

6.7%  93.3%  21.2% 

3.5% 

24.7% 

0 

12.6% 

B 

6.98E-04  0.0%  100.0%  42.0% 

40.5% 

-3.5% 

45.7% 

0 

40.5% 

C 

5.93E-04  82.2%  17.8%  8.7% 

7.1% 

-18.0% 

22.4% 

0 

7.1% 

AVE 

7.23E-04  52.6%  47.4%  21.0% 

20.1% 

-6.0% 

30.9% 

0 

20.1% 

0.28 

A 

7.35E-03  68.4%  32.6%  5.6% 

5.6% 

39.6%  60.4%  3.7% 

-0.3% 

6.7% 

0 

5.6% 

B 

7.39E-03  0.0%  100.0%  4.7% 

4.7% 

0.3% 

6.0% 

0 

4.7% 

C 

7.51E-03  30.4%  69.6%  1.4% 

1.4% 

1.9% 

3.9% 

0 

1.4% 

AVE 

7.37E-03  32.9%  67.4%  3.9% 

3.9% 

0.6% 

5.5% 

0 

3.9% 

0.41 

A 

1.73E-02  56.3%  44.7%  2.5% 

2.4% 

21.4%  78.6%  1.9% 

-4.2% 

3.1% 

2.9% 

3.7% 

B 

1.88E-02  0.0%  100.0%  2.9% 

3.1% 

4.2% 

3.6% 

2.2% 

3.7% 

C 

1.85E-02  42.8%  57.2%  0.6% 

0.6% 

2.7% 

2.0% 

1.8% 

1.9% 

AVE 

1.80E-02  33.0%  67.3%  2.0% 

2.0% 

0.9% 

2.9% 

2.3% 

3.1% 

Note: 

(1)  DRE  for  sinkage  in  ONR  paper  figures  for  three  institutes  is 

2(AFP  +  AAP)  .  ,  „  2  2(AFP  +  AAP) 

cr  = -  instead  ot  cr  =  — - - 

L  Fr 2  2 L 

(2)  Uncertainties  in  terms  of  percentages  remain  the  same  for  both  above  DREs, 
while  uncertainty  magnitudes  are  different  as  long  as  we  assume  that  Fr  number 
has  no  uncertainty  for  simplicity. 

(3)  X1  for  A,  B,  and  C  from  Interpolated  data  files  (dtmbintp.tec,  seanintp.tec, 
iihrintp.tec) 

(4)  All  B\  and  Pi  from  [16] 

(5) B2X  and  Pi  for  DTMB  (Fr=0.28  and  0.41)  are  best  estimates  from  [16] 

(6)  All  U x  based  on  mean  values  from  Interpolated  data  files  (dtmb  intp.tec, 
sean  intp.tec,  iihr  intp.tec) 
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Table  C2.  Facility  certification  or  biases  for  trim  results  z 


Fr  & 

Facility 

N-order  level  (%  Xt) 

MxN-order  level  (%  X  ) 

Facility  Certification 

or  biases  (%  X  ) 

Pi  UXt 

Ux, 

4  4  ux 

A 

uDl 

UFB, 

uTj 

0.10 

A 

-2.05E-04  64.5%  35.5%  14.4% 

7.2% 

1.4%  98.6%  24.3% 

-49.9% 

25.3% 

43.0% 

43.6% 

B 

-6.13E-04  0.0%  100.0%  32.0% 

48.0% 

49.9% 

53.8% 

0 

48.0% 

C 

-6.22E-04  50.8%  49.2%  10.2% 

15.6% 

52.3% 

28.8% 

43.7% 

46.4% 

AVE 

-4.09E-04  38.4%  61.6%  18.9% 

23.6% 

17.4% 

36.0% 

28.9% 

46.0% 

0.28 

A 

-3.90E-03  54.7%  46.3%  2.8% 

2.8% 

15.0%  85.0%  2.7% 

1.7% 

3.9% 

0 

2.8% 

B 

-3.77E-03  0.0%  100.0%  4.7% 

4.6% 

-1.7% 

5.4% 

0 

4.6% 

C 

-4.75E-03  36.1%  63.9%  1.8% 

2.3% 

23.9% 

3.5% 

23.6% 

23.7% 

AVE 

-3.83E-03  30.3%  70.1%  3.1% 

3.2% 

8.0% 

4.3% 

7.9% 

10.4% 

0.41 

A 

1.36E-02  38.1%  61.9%  1.5% 

1.5% 

28.0%  72.0%  0.9% 

-1.6% 

1.7% 

0 

1.5% 

B 

1.40E-02  0.0%  100.0%  0.9% 

0.9% 

1.6% 

1.2% 

1.0% 

1.3% 

C 

1.06E-02  4.1%  95.9%  1.8% 

1.4% 

-23.1% 

1.6% 

23.0% 

23.0% 

AVE 

1.38E-02  14.1%  85.9%  1.4% 

1.2% 

-7.7% 

1.5% 

8.0% 

8.6% 

Note: 

(1)  DRE  for  trim  in  ONR  paper  figures  for  three  institutes  is 

t_  2(AAP  -  AFP)  Mo(r  =  ^(AAP-AFP} 

L  Fr2  L 

(2)  Uncertainties  in  terms  of  percentages  remain  the  same  for  both  above  DREs, 
while  uncertainty  magnitudes  are  different  as  long  as  we  assume  that  Fr  number 
has  no  uncertainty  for  simplicity. 

(3)  X1  for  A,  B,  and  C  from  Interpolated  data  files  (dtmbintp.tec,  seanintp.tec, 
iihrintp.tec) 

(4)  All  B\  and  P2  from  [16] 

(5 ) B\  and  P2  for  DTMB  (Fr=0.28)  are  best  estimates  from  [16] 

(6)  All  U x  based  on  mean  values  from  Interpolated  data  files  (dtmb  intp.tec, 
sean  intp.tec,  iihr  intp.tec) 
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Table  C3.  Facility  certification  or  biases  for  wave  profile  results  C, 


Fr  & 

Facility 

N-order  level  (%  Dx  ) 

MxN-order  level  (%  Dx  ) 

Facility  Certification 

or  biases  (%  Dx  ) 

Xt  DXi  B%  Pi  UXi 

Ux, 

B-  P-  U- 

X  1  X  X 

A 

uDi 

uFBi 

UTj 

0.28 

A 

1.82E-02  1.89E-02  64.5%  35.5%  3.5% 

3.2% 

87.0%  13.0%  2.7% 

-4.4% 

4.2% 

1.4% 

3.5% 

B 

2.00E-02  2.09E-02  100.0%  0.0%  4.2% 

4.2% 

4.4% 

5.0% 

0 

4.2% 

C 

2.03E-02  2.22E-02  83.7%  16.3%  3.4% 

3.7% 

5.9% 

4.5% 

3.8% 

5.3% 

AVE 

1.91E-02  2.07E-02  82.7%  17.3%  3.7% 

3.7% 

2.0% 

4.6% 

1.7% 

4.3% 

0.41 

A 

3.02E-02  3.62E-02  64.5%  35.5%  1.8% 

2.0% 

81.5%  18.5%  1.4% 

3.1% 

2.4% 

1.9% 

2.8% 

B 

2.82E-02  2.46E-02  100.0%  0.0%  2.6% 

1.9% 

-3.1% 

2.4% 

2.0% 

2.8% 

C 

3.03E-02  3.86E-02  81.6%  18.4%  2.0% 

2.3% 

3.3% 

2.7% 

1.9% 

3.0% 

AVE 

2.92E-02  3.32E-02  82.0%  18.0%  2.1% 

2.1% 

1.1% 

2.5% 

1.9% 

2.9% 

X,  is  the  maximum  elevation  on  the  wave  profile. 


Note: 

(1)  Xt  from  Data  files  before  interpolation  (dtmb_dpwpro28.tec, 
sean_dpwpro28.tec,  iihr_dpwpro28.tec) 

(2)  Dx  for  A  from  DTMB’s  website:  http://www50.dt.navy.mil/5415/profile.html, 
(0.15  inch  for  all  points) 

(3)  Dx  for  B  from  IIHR  Report  421  (p.53)  for  INSEAN 

(4)  Dx  for  C  from  IIHR’s  raw  data  files  (wpro2801.tec,  wpro2802.tec, 

wpro2803.tec,  wpro4101.tec,  wpro4102.tec,  wpro4103.tec) 

(5)  All  B2x  and  PXj  from  [16] 

(6)  All  JJ x  from  [16] 
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Table  C4.  Facility  certification  or  biases  for  wave  elevation  results  C,  at  cut  y=0.324 


Fr  & 

Facility 

N-order  level  (%  Dx  ) 

MxN-order  level  (%  Dx  ) 

Facility  Certification 

or  biases  (%  Dx  ) 

Xt  DXi  B\  Pi  UXt 

uXi 

4  4  u-x 

A 

uDi 

UFB, 

uTj 

0.28 

4.19E-03  1.17E-02  76.6%  23.4%  2.7% 

2.8% 

71.9%  28.1%  1.8% 

-4.5% 

3.4% 

3.0% 

4.1% 

B 

5.20E-03  1.09E-02  64.9%  35.1%  2.4% 

2.3% 

4.5% 

3.0% 

3.4% 

4.1% 

3.79E-03  1.11E-02  59.0%  41.0%  3.4% 

3.4% 

-8.1% 

3.8% 

7.1% 

7.9% 

4.70E-03  1.12E-02  66.8%  33.2%  2.9% 

2.9% 

-2.7% 

3.4% 

4.5% 

5.4% 

Xt  is  the  maximum  elevation  on  the  cut  y=0.324. 


Note: 

(1)  Xt  from  Interpolated  data  files  (dtmb_int324.tec,  insean_int324.tec, 
iihr_int324.tec) 

(2)  Dx  for  A  from  Interpolated  data  files  (dtmb  intpat.tec) 

(3)  Dx  for  B  from  IIHR  Report  421  (p.54)  for  INSEAN 

(4)  Dx  for  C  from  (Joe  Longo)  IIHR’s  raw  data  files  (a0_101.dat  zone  =  “Steady”) 

(5)  All  B\  and  Pi  from  [16] 

(6)  All  UXj  from  [16] 

(7 ) BX  and  Px  for  DTMB  (Fr=0.28)  are  best  estimates  from  [16] 
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Table  C5.  Facility  certification  or  biases  for  wave  elevation  results  C, 


Fr& 

Facility 

N-order  level  (%  D x  ) 

MxN-order  level  (%  Dx  ) 

Facility  Certification 

or  biases  (%  Dx  ) 

Ar, 

4,.  Pi  UXi 

uXl 

4  4  ux 

A 

uDi 

uFBi 

UTj 

0.28 

A 

1.17E-02 

76.6%  23.4%  2.7% 

2.8% 

71.9%  28.1%  1.8% 

3.4%  3.4%  0  2.8% 

B 

1.09E-02 

64.9%  35.1%  2.4% 

2.3% 

3.4%  3.0%  1.6%  2.8% 

C 

1.11E-02 

59.0%  41.0%  3.4% 

3.4% 

2.9%  3.8%  0  3.4% 

AVE 

1.12E-02 

66.8%  33.2%  2.9% 

2.9% 

3.2%  3.4%  0.6%  3.0% 

Note: 

(1)  Dx  for  A  from  Interpolated  data  files  (dtmbintpat.tec) 

(2)  Dx  for  B  from  IIHR  Report  421  (p.54)  for  INSEAN 

(3)  Dx  for  C  from  (Joe  Longo)  IIHR’s  raw  data  files  (a0_101.dat  zone  =  “Steady”) 

(4)  All  B2Xj  and  PXj  from  [16] 

(5)  All  UXj  from  [16] 

(6 ) BX  and  Px  for  DTMB  (Fr=0.28)  are  best  estimates  from  [16] 

(7) A11  Di  calculated  from  Interpolated  data  files  (dtmb  intpat.tec,  insean  intpat.tec, 
iihrintpat.tec) 
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Table  C6.  Facility  certification  or  biases  for  nominal  wake  velocity  results  at  Fr=0.28 


N-order  level  (%  Dx  ) 

Facility  Certification 

Fr  & 

MxN-order  level  (%  Dx  ) 

or 

biases  (%  D x  ) 

Facility 

LTR 

HTR 

SQRT  (LTRA2+HTRA2)/2 

DXl 

Br  P2 

Xi  Xi 

uXt 

B 2  P2 

Xi  Xi 

uXl 

K 

P 2 

Xi 

uXi 

uXi 

B 2  P- 

X  1  X 

Ux 

A 

U D, 

uFBi 

UTj 

A 

1.00E+00 

65.4%  34.4% 

1.6% 

74.5%  26.5% 

12.5% 

73.6% 

26.4% 

6.3% 

6.3% 

1.0% 

7.0% 

0 

6.3% 

0.28 

B 

1.02E+00 

47.8%  52.2% 

0.4% 

60.4%  39.6% 

0.4% 

53.5% 

46.5% 

0.3% 

0.3% 

73.6%  26.4% 

3.1% 

1.0% 

3.2% 

0 

0.3% 

U 

C 

1.01E+00 

99.8%  0.2% 

1.2% 

99.2%  0.8% 

3.1% 

99.3% 

0.7% 

1.7% 

1.7% 

2.5% 

3.6% 

0 

1.7% 

AVE 

1.01E+00 

71.0%  28.9% 

1.1% 

78.0%  22.3% 

5.3% 

75.5% 

24.5% 

2.8% 

2.7% 

1.5% 

4.6% 

0 

2.7% 

A 

9.62E-02 

54.3%  46.7% 

2.9% 

43.5%  56.5% 

6.5% 

45.2% 

54.8% 

3.6% 

3.0% 

3.6% 

3.4% 

1.2% 

3.2% 

0.28 

B 

1.19E-01 

21.2%  78.8% 

1.9% 

15.9%  84.1% 

2.7% 

17.6% 

82.4% 

1.6% 

1.7% 

38.4%  61.6% 

1.7% 

3.6% 

2.4% 

2.7% 

3.2% 

V 

C 

1.33E-01 

99.7%  0.3% 

5.5% 

93.3%  6.7% 

3.8% 

97.7% 

2.3% 

3.3% 

3.8% 

6.0% 

4.2% 

4.3% 

5.8% 

AVE 

1.16E-01 

58.4%  41.9% 

3.4% 

50.9%  49.1% 

4.3% 

53.5% 

46.5% 

2.9% 

2.8% 

4.4% 

3.3% 

2.7% 

4.0% 

A 

1.47E-01 

65.3%  34.7% 

6.5% 

44.6%  35.4% 

3.7% 

63.3% 

36.7% 

3.6% 

3.3% 

8.6% 

3.7% 

7.8% 

8.5% 

0.28 

B 

1.65E-01 

79.1%  20.9% 

1.0% 

87.9%  12.1% 

0.9% 

83.1% 

16.9% 

0.6% 

0.6% 

64.1%  35.9% 

1.7% 

8.6% 

1.8% 

8.4% 

8.5% 

W 

C 

1.79E-01 

99.9%  0.1% 

4.1% 

99.2%  0.8% 

4.5% 

99.5% 

0.5% 

3.0% 

3.3% 

8.8% 

3.7% 

7.9% 

8.6% 

AVE 

1.64E-01 

81.4%  18.6% 

3.8% 

77.2%  16.1% 

3.0% 

82.0% 

18.0% 

2.4% 

2.4% 

8.7% 

3.1% 

8.1% 

8.5% 

Note: 

(1)  Dx  for  A  from  DTMB  Raw  data  files:  (dtmb  dpnwake.tec) 

(2)  Dx  for  B  from  IIHR  Report  421  (p.56)  for  INSEAN,  but  dynamic  range  for  U  is  divided  by  2. 

(3)  Dx  for  C  from  IIHR  Raw  data  document 
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(4)  All  B\  and  P],  from  [16] 

(5)  All  U y_  from  [16] 

(6 ) B\  and  P Y2  for  DTMB  are  best  estimates  from  [16] 

(7)  All  Z>  calculated  from  Interpolated  data  files  (dtmbdpINTnwake.tec,  sean  dpINTnwake.tec, 
iihrdpINT  nwake.tec) 
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Figure  Cl  Residuary  resistance  results  and  uncertainty  bands  at  Fr=0.1,  0.28  and 

0.41 


Note: 

(1)  DTMB  data  file:  \23rdONRdata\Resistance\dtmb_dpcr.tec 

(2)  INSEAN  data  file:  \23rdONRdata\Resistance\sean_dpcr.tec 

(3)  IIHR  data  file:  \23rdONRdata\Resistance\iihr_dpcr.tec 

(4)  Mean  data  file:  \23rdONRdata\Resistance\meanAB.txt 

(5)  DTMB  uncertainty  file:  \23rdONRdata\Resistance\dtmb_dpcr_un.dat 

(6)  INSEAN  uncertainty  file:  \23rdONRdata\Resistance\sean_dpcr_un.dat 

(7)  IIHR  uncertainty  file:  \23rdONRdata\Resistance\iihr_dpcr_un.dat 
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Figure  C2  Sinkage  results  and  uncertainty  bands  at  Fr=0.1,  0.28  and  0.41 


Note: 

(1)  DTMB  data  file:  \23rdONRdata\Sinkage_Trim\dtmb_dpst.tec 

(2)  INSEAN  data  file:  \23rdONRdata\Sinkage_Trim\sean_dpst.tec 

(3)  IIHR  data  file:  \23rdONRdata\Sinkage_Trim\iihr_dpst.tec 

(4)  Mean  data  file:  \23rdONRdata\Sinkage_Trim\meanAB.txt 

(5)  DTMB  uncertainty  file:  \23rdONRdata\Sinkage_Trim\dtmb_dpcr_sink_un.dat 

(6)  INSEAN  uncertainty  file:  \23rdONRdata\Sinkage_Trim\sean_dpcr_sink_un.dat 

(7)  IIHR  uncertainty  file:  \23rdONRdata\Sinkage_Trim\iihr_dpcr_sink_un.dat 
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Figure  C3  Trim  results  and  uncertainty  bands  at  Fr=0.1,  0.28  and  0.41 


Note: 

(1)  DTMB  data  file:  \23rdONRdata\Sinkage_Trim\dtmb_dpst.tec 

(2)  INSEAN  data  file:  \23rdONRdata\Sinkage_Trim\sean_dpst.tec 

(3)  IIHR  data  file:  \23rdONRdata\Sinkage_Trim\iihr_dpst.tec 

(4)  Mean  data  file:  \23rdONRdata\Sinkage_Trim\meanAB.txt 

(5)  DTMB  uncertainty  file:  \23rdONRdata\Sinkage_Trim\dtmb_dpcr_trim_un.dat 

(6)  INSEAN  uncertainty  file:  \23rdONRdata\Sinkage_Trim\sean_dpcr_trim_un.dat 

(7)  IIHR  uncertainty  file:  \23rdONRdata\Sinkage_Trim\iihr_dpcr_trim_un.dat 
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Figure  C4  Wave  profile  results  and  uncertainty  bands  at  ^max  (Fr=0.28) 


Note: 

(1)  DTMB  data  file:  \23rdONRdata\Wave_Profiles\dtmb_dpwpro28.tec 

(2)  INSEAN  data  file:  \23rdONRdata\Wave_Profiles\sean_dpswpro28.tec 

(3)  IIHR  data  file:  \23rdONRdata\Wave_Profiles\iihr_dpwpro28.tec 

(4)  Mean  data  file:  \23rdONRdata\Wave_Profiles\mean28AB.txt 

(5)  DTMB  uncertainty  file:  \23rdONRdata\Wave_Profiles\dtmb_dp28_un.dat 

(6)  INSEAN  uncertainty  file:  \23rdONRdata\Wave_Profiles\sean_dp28_un.dat 

(7)  IIHR  uncertainty  file:  \23rdONRdata\Wave_Profiles\iihr_dp28_un.dat 


56 


Figure  C5  Wave  profile  results  and  uncertainty  bands  at  (Fr=0.41) 


Note: 

(1)  DTMB  data  file:  \23rdONRdata\Wave_Profiles\dtmb_dpwpro41.tec 

(2)  INSEAN  data  file:  \23rdONRdata\Wave_Profiles\sean_dpswpro41.tec 

(3)  IIHR  data  file:  \23rdONRdata\Wave_Profiles\iihr_dpwpro41.tec 

(4)  Mean  data  file:  \23rdONRdata\Wave_Profiles\mean41AB.txt 

(5)  DTMB  uncertainty  file:  \23rdONRdata\Wave_Profiles\dtmb_dp41_un.dat 

(6)  INSEAN  uncertainty  file:  \23rdONRdata\Wave_Profiles\sean_dp41_un.dat 

(7)  IIHR  uncertainty  file:  \23rdONRdata\Wave_Profiles\iihr_dp41_un.dat 
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Figure  C6  Wave  elevation  results  and  uncertainty  bands  for  cut  y=0.324  (Fr=0.28) 

Note: 

(1)  DTMB  data  file:  \23rdONRdata\Wave_Elevations\dtmb_int324.tec 

(2)  INSEAN  data  file:  \23rdONRdata\Wave_Elevations\sean_int324.tec 

(3)  IIHR  data  file:  \23rdONRdata\Wave_Elevations\iihr_int324.tec 

(4)  Mean  data  file:  \23rdONRdata\Wave_Elevations\meanAB.txt 

(5)  DTMB  uncertainty  file:  \23rdONRdata\Wave_Elevations\dtmb_dp28_un.dat 

(6)  INSEAN  uncertainty  file:  \23rdONRdata\Wave_Elevations\sean_dp28_un.dat 

(7)  IIHR  uncertainty  file:  \23rdONRdata\Wave_Elevations\iihr_dp28_un.dat 


58 


DTMB  541 5  Far  Field  Wave  Elevation  Results 


INSEAN  2340A  Far  Field  Wave  Elevation  Results 


IIHR  5512  Far  Field  Wave  Elevation  Results 


Mean  Far  Field  Wave  Elevation  Results 


Figure  C7  Wave  elevation  results  for  far  field  and  the  mean  result  (Fr=0.28) 

Note: 

(1)  DTMB  data  file:  \23rdONRdata\Wave_Elevations\dtmb_intpat.tec 

(2)  INSEAN  data  file:  \23rdONRdata\Wave_Elevations\insean_intpat.tec 

(3)  IIHR  data  file:  \23rdONRdata\Wave_Elevations\iihr_intpat.tec 

(4)  Mean  data  file:  \23rdONRdata\Wave_Elevations\mean_Zeta_AB.dat 
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D.  for  INSEAN  2340A  Far  Field  Wave  Elevation  Results 
Dj  =  3.4%,  UDi  =  3.0%  (  %D. ) 


D.  for  IIHR  5512  Far  Field  Wave  Elevation  Results 
Di  =  2.9%,  UDi  =  3.8%  (%DC) 


Figure  C8  Contours  of  difference  Di  of  wave  elevation  results  (Fr=0.28) 


Note: 

(1)  DTMB  data  file:  \23rdONRdata\Wave_Elevations\DTMB_DZeta_AB.tec 

(2)  INSEAN  data  file:  \23rdONRdata\Wave_Elevations\SEAN  DZeta  AB.tec 

(3)  IIHR  data  file:  \23rdONRdata\Wave_Elevations\IIHR_DZeta_AB.tec 
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Certification  of  Far  Field  Wave  Elevation  Results  (DTMB) 
Mean  =  (A+B)/2 


UDi  =  3.4%D 


Certification  of  Far  Field  Wave  Elevation  Results  (INSEAN) 
Mean  =  (A+B)/2 


61 


Figure  C9  3D  surface  of  difference  Di  of  wave  elevation  results  (Fr=0.28) 

Note: 

(1)  DTMB  data  file:  \23rdONRdata\Wave_Elevations\DTMB_DZeta_AB.tec 

(2)  INSEAN  data  fde:  \23rdONRdata\Wave_Elevations\SEAN  DZeta  AB.tec 

(3)  IIHR  data  file:  \23rdONRdata\Wave_Elevations\IIHR_DZeta_AB.tec 
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Figure  CIO  Nominal  wake  velocity  results  and  the  mean  velocity  result  (Fr=0.28) 
Note: 

(1)  DTMB  data  fde:  \23rdONRdata\Nominal  Wake\DTMB  dpINTnwake.dat 

(2)  INSEAN  data  fde:  \23rdONRdata\Nominal  Wake\SEAN_dpINTnwake.dat 

(3)  IIHR  data  file:  \23rdONRdata\Nominal  WakeUIHR dpINTnwake.dat 

(4)  Mean  data  file:  \23rdONRdata\Nominal_Wake\MeanAB_UVW.dat 
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Figure  Cll  Contours  of  difference  Di  of  nominal  wake  velocity  results  (Fr=0.28) 
Note: 

(1)  DTMB  data  file:  \23rdONRdata\Nominal_Wake\MeanAB_DA_perc.dat 

(2)  INSEAN  data  file:  \23rdONRdata\Nominal_Wake\MeanAB_DB_perc.dat 

(3)  IIHR  data  file:  \23rdONRdata\Nominal_Wake\MeanAB_DC_perc.dat 
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Certification  of  nominal  wake  velocity  component  W  (IIHR) 
Mean  =  (A+B)/2 


Figure  C12  3D  surface  of  difference  Di  of  nominal  wake  velocity  results  (Fr=0.28) 
Note: 

(1)  DTMB  data  file:  \23rdONRdata\Nominal  Wake\MeanAB  DA.dat 

(2)  INSEAN  data  fde:  \23rdONRdata\Nominal  Wake\MeanAB  DB.dat 

(3)  IIHR  data  file:  \23rdONRdata\Nominal  Wake\MeanABC_D.dat 
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